An analysis of overlapping Schwarz method for a weakly coupled system of singularly perturbed convection-diffusion equations

In this article, we have developed an overlapping Schwarz method for a weakly coupled system of convection-diffusion equations. The method splits the original domain into two overlapping subdomains. A hybrid difference scheme is proposed in which on the boundary layer region, we use the central finite difference scheme on a uniform mesh, whereas on the nonlayer region, we use the mid-point difference scheme on a uniform mesh. It is shown that the numerical approximations converge in the maximum norm to the exact solution. We have proved that, when appropriate subdomains are used, the method produces almost second-order convergence. Furthermore, it is shown that two iterations are sufficient to achieve the expected accuracy. Numerical examples are presented to support the theoretical results. The main advantage of this method used with the proposed scheme is that it reduces iteration counts very much and easily identifies in which iteration the Schwarz iterate terminates.     

A control volume finite element method for three-dimensional three-phase flows

A novel control volume finite element method with adaptive anisotropic unstructured meshes is presented for three-dimensional three-phase flows with interfacial tension. The numerical framework consists of a mixed control volume and finite element formulation with a new P₁DG-P₂ elements (linear discontinuous velocity between elements and quadratic continuous pressure between elements). A “volume of fluid” type method is used for the interface capturing, which is based on compressive control volume advection and second-order finite element methods. A force-balanced continuum surface force model is employed for the interfacial tension on unstructured meshes. The interfacial tension coefficient decomposition method is also used to deal with interfacial tension pairings between different phases. Numerical examples of benchmark tests and the dynamics of three-dimensional three-phase rising bubble, and droplet impact are presented. The results are compared with the analytical solutions and previously published experimental data, demonstrating the capability of the present method.     

Stochastic sensitivity analysis of coupled acoustic-structural systems under non-stationary random excitations

This paper presents a new sensitivity analysis method for coupled acoustic–structural systems subjected to non-stationary random excitations. The integral of the response power spectrum density (PSD) of the coupled system is taken as the objective function. The thickness of each structural element is used as a design variable. A time-domain algorithm integrating the pseudo excitation method (PEM), direct differentiation method (DDM) and high precision direct (HPD) integration method is proposed for the sensitivity analysis of the objective function with respect to design variables. Firstly, the PEM is adopted to transform the sensitivity analysis under non-stationary random excitations into the sensitivity analysis under pseudo transient excitations. Then, the sensitivity analysis equation of the coupled system under pseudo transient excitations is derived based on the DDM. Moreover, the HPD integration method is used to efficiently solve the sensitivity analysis equation under pseudo transient excitations in a reduced-order modal space. Numerical examples are presented to demonstrate the validity of the proposed method.     

A proper orthogonal decomposition variational multiscale meshless interpolating element-free Galerkin method for incompressible magnetohydrodynamics flow

In the recent decade, the meshless methods have been handled for solving most of PDEs due to easiness of the meshless methods. One of the popular meshless methods is the element-free Galerkin (EFG) method that was first proposed for solving some problems in the solid mechanics. The test and trial functions of the EFG are based on the special basis. Recently, some modifications have been developed to improve the EFG method. One of these improvements is the variational multiscale EFG procedure. In the current article, the shape functions of interpolation moving least squares approximation have been applied to the variational multiscale EFG technique for solving the incompressible magnetohydrodynamics flow. In order to reduce the elapsed CPU time of simulation, we employ a reduced-order model based on the proper orthogonal decomposition technique. The current combination can be referred to as the reduced-order variational multiscale EFG technique. To illustrate the reduction in CPU time used as well as the efficiency of the proposed method, we applied it for the two-dimensional cases.     

Momentum and heat transfer of a special case of the unsteady stagnation-point flow

This paper investigates the unsteady stagnation-point flow and heat transfer over a moving plate with mass transfer,which is also an exact solution to the unsteady Navier-Stokes(NS)equations.The boundary layer energy equation is solved with the closed form solutions for prescribed wall temperature and prescribed wall heat flux conditions.The wall temperature and heat flux have power dependence on both time and spatial distance.The solution domain,the velocity distribution,the flow field,and the temperature distribution in the fluids are studied for different controlling parameters.These parameters include the Prandtl number,the mass transfer parameter at the wall,the wall moving parameter,the time power index,and the spatial power index.It is found that two solution branches exist for certain combinations of the controlling parameters for the flow and heat transfer problems.The heat transfer solutions are given by the confluent hypergeometric function of the first kind,which can be simplified into the incomplete gamma functions for special conditions.The wall heat flux and temperature profiles show very complicated variation behaviors.The wall heat flux can have multiple poles under certain given controlling parameters,and the temperature can have significant oscillations with overshoot and negative values in the boundary layers.The relationship between the number of poles in the wall heat flux and the number of zero-crossing points is identified.The difference in the results of the prescribed wall temperature case and the prescribed wall heat flux case is analyzed.Results given in this paper provide a rare closed form analytical solution to the entire unsteady NS equations,which can be used as a benchmark problem for numerical code validation.     

A new database and benchmark of the bond energies of noble-gas-containing molecules

We have developed a new database of structures and bond energies of 59 noble-gas-containing molecules. The structures were calculated by CCSD(T)/aug-cc-pVTZ methods and the bond energies were obtained using the CCSD(T)/complete basis set method. Many wavefunction-based and density functional theory methods have been benchmarked against the 59 accurate bond energies. Our results show that the MPW1B95, B2GP-PLYP, and DSD-BLYP functionals with the aug-cc-pVTZ basis set excel in predicting the bond energies of noble-gas molecules with mean unsigned errors (MUEs) of 2.0 to 2.1 kcal/mol. When combinations of Dunning's basis sets are used, the MPW1B95, B2GP-PLYP, DSD-BLYP, and BMK functionals give significantly lower MUEs of 1.6 to 1.9 kcal/mol. Doubly hybrid methods using B2GP-PLYP and DSD-BLYP functionals and MP2 calculation also provide satisfactory accuracy with MUEs of 1.4 to 1.5 kcal/mol. If the Ng bond energies and the total atomization energies of a group of 109 main-group molecules are considered at the same time, the MPW1B95/aug-cc-pVTZ single-level method (MUE = 2.7 kcal/mol) and the B2GP-PLYP and DSD-PLYP functionals with combinations of basis sets or using the doubly hybrid method (MUEs = 1.9-2.2 kcal/mol) give the overall best result.     

Multi-variable special polynomials using an operator ordering method

Using an operator ordering method for some commutative superposition operators, we introduce two new multi-variable special polynomials and their generating functions, and present some new operator identities and integral formulas involving the two special polynomials. Instead of calculating complicated partial differential, we use the special polynomials and their generating functions to concisely address the normalization, photocount distributions and Wigner distributions of several quantum states that can be realized physically, the results of which provide real convenience for further investigating the properties and applications of these states.     

A new non-perturbative approach in quantum mechanics for time-independent Schr?dinger equations

A new non-perturbative approach is proposed to solve time-independent Schr?dinger equations in quantum mechanics.It is based on the homotopy analysis method(HAM)that was developed by the author in 1992 for highly nonlinear equations and has been widely applied in many fields.Unlike perturbative methods,this HAM-based approach has nothing to do with small/large physical parameters.Besides,convergent series solution can be obtained even if the disturbance is far from the known status.A nonlinear harmonic oscillator is used as an example to illustrate the validity of this approach for disturbances that might be one thousand times larger than the possible superior limit of the perturbative approach.This HAM-based approach could provide us rigorous theoretical results in quantum mechanics,which can be directly compared with experimental data.Obviously,this is of great benefit not only for improving the accuracy of experimental measurements but also for validating physical theories.